Kapellasite: a kagome quantum spin liquid with competing interactions 
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Magnetic susceptibility, NMR, /iSR, and inelastic neutron scattering measurements show that 
kapellasite, Cu3Zn(OH)6Cl2, a geometrically frustrated spin-1/2 kagome antiferromagnet polymor- 
phic with herbertsmithite, is a gapless spin liquid showing unusual dynamic short-range correlations 
of non-coplanar cuboc2 type which persist down to 20 mK. The Hamiltonian is determined from a 
fit of a high-temperature series expansion to bulk susceptibility data and possesses competing ex- 
change interactions. The magnetic specific heat calculated from these exchange couplings is in good 
agreement with experiment. The temperature dependence of the magnetic structure factor and the 
muon relaxation rate are calculated in a Schwinger-boson approach and compared to experimental 
results. 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.-s, 78.70.Nx 



The quantum spin liquid (QSL) is one of the most elu- 
sive states in condensed matter physics, even after two 
decades of intense studies [1| , with implications ranging 
from geometrically or exchange frustrated magnetic sys- 
tems to unconventional superconductivity. A promising 
model system for the observation of the QSL state is the 
5"= 1/2 kagome antiferromagnet (KAFM), which is char- 
acterized by quantum spins on a geometrically frustrated 
lattice of corner-sharing equilateral triangles. The first 
suitable model material, the mineral herbertsmithite, is 
classified as a gapless spin liquid 043' ^^"^ still lacks a 
complete understanding. Its discovery triggered a burst 
of experimental and theoretical investigations in the field 
of quantum spin liquids in two dimensions and challenges 
the most advanced theoretical approaches |9l4l2|. 

A new experimental realization of the 5"= 1/2 KAFM 
model was recently synthesized [l3|, |lj|, the mineral 
kapellasite, Cu3Zn(OH)6Cl2, a polymorph of herbert- 
smithite, i.e., they have the same chemical formula but 
different crystallographic structures. The layered struc- 
ture of kapellasite prevents the Cu/Zn disorder present in 
these materials from either coupling the kagome planes or 
creating free spins, in contrast to the more 3-dimcnsional 
herbertsmithite. In this Letter, we use inelastic neu- 
tron scattering and muon spin relaxation (^uSR) down 
to 20 mK to show that kapellasite is a gapless spin liquid 
without any sign of long-range magnetic order or spin 
freezing. In contrast to the featureless excitation spec- 
trum of herbertsmithite, kapellasite exhibits novel dy- 
namic short-range correlations with a well-defined wave- 
vector dependence (see Fig. [T]), reminiscent of the non- 



coplanar twelve-sublattice cuboc2 magnetic structure 
[l5|, where, due to a particular set of fcrro- and anti- 
ferromagnetic competing interactions, spins on a given 
hexagon are coplanar with neighboring spins forming an 
angle of 60° with each other; the third-neighbor spins 
across hexagons are antiparallel. 

Kapellasite crystallizes in the PimX trigonal space 
group (#164) with lattice parameters a=6.33 and 




FIG. 1. (Color online) Neutron scattering intensity, S{Q,E), 
on a logarithmic scale as a function of wave vector Q and 
energy E of kapellasite at a temperature of 0.1 K measured 
with an incoming neutron energy of 3.27 meV. The diffuse 
inelastic magnetic scattering is centered at (3=0.5 and 2.0 
A~^. On the right are shown the spin directions (bottom) on 
the kagome lattice (middle) of the cuboc2 magnetic structure. 




cuboci 



VsxVs' 






cuboc2 


* 


_ ferro 


k 


1 1 


k=0 



-2-1012 

Jo 



FIG. 2. (Color online) (a) Kagome plane of kapellasite with 
Cu^^ S — 1/2 spins (blue), non- magnetic Zn^"'" ion (green), 
and exchange interactions (red), (b) Classical J2-Jd phase 
diagram for ferromagnetic nearest-neighbor coupling, Ji — 
— 1, and Ja = 0. The best two-parameter HT-series fits to 
susceptibility data are shown with red points. 



0=5.70 A. The S—1/2 kagome lattice of undistorted 
corner-sharing triangles is obtained by regularly dop- 
ing a two-dimensional (2D) triangular Cu^+ metal-site 
sublattice with diamagnetic Zn-^+ ions [see Fig. [DJa)]. 
Inductively coupled plasma mass spectrometry showed 
that our sample has not the ideal Cu/Zn 3:1 stoi- 
chiometry and '^^Cl NMR showed some Cu/Zn inter- 
site mixing [16|. The actual chemical formula consis- 
tently determined with neutron powder diffraction |14{ 
is (Cuo.73Zno.27)3(Zno.88Cuo.i2)(OH)6Cl2, with 27% Zn 
on the Cu-sites of the kagome lattice and 12% Cu on the 
"hexagonal" Zn site. The Cu/Zn mixing leads to some 
disorder within the kagome planes, but cannot induce 
any coupling between the planes, which occurs only via 
very weak 0-H-Cl hydrogen bonds [l3j. Kapellasite is 



therefore remarkably two-dimensional. 

The DC magnetic susceptibility measured in a SQUID 
in a magnetic field of 5 T exhibits a Curie- Weiss behavior 
Xuc{T) = C/{T — 6cw) above ^^50 K with a weak ferro- 
magnetic effective coupling, 9cw = 9.5± 1 K [Fig.[31Ja)]. 
For r < 50 K, the susceptibility becomes Curie-like, with 
neither any sign of a magnetic transition nor a FC/ZFC 
opening. A similar Curie-tail is observed in many other 
compounds and is often attributed to a defect contribu- 
tion which exceeds the intrinsic susceptibility la,[l4|- The 
use of a local probe such as NMR is decisive in this con- 
text: We find that the local susceptibility xioc extracted 
from the ^^Cl NMR shift of an oriented powder sample 
with magnetic fields up to 13 T applied perpendicular to 
the kagome planes, scales with xdc over the whole tem- 
perature range 2-287 K. This clearly shows that isolated 
spin defects are not present in kapellasite in contrast with 
the behavior seen in herbertsmithite, where Cu/Zn mix- 
ing creates paramagnetic-like spins. Thus, the intrinsic 
susceptibility of kapellasite displays genuine divergent be- 
havior pointing to a gapless spin state, rather than the 
downturn seen in other KAFM compounds. 
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FIG. 3. (Color online) (a) Local xioc and macroscopic xuc 
magnetic susceptibilities from NMR and DC SQUID mea- 
surements, respectively. Inset: Fit of high-temperature series 
expansion (red line) down to 20 K. (b) Total specific heat per 
spin C™^ measured in zero field (symbols), magnetic specific 
heat calculated with the J\~Jd model (red solid line) rescaled 
by a factor 0.88 to account for mass uncertainties and missing 
entropy. The phonon contribution estimated from a high-T 
fit (purple dashed line) is negligible below 10 K. 



A reasonable first theoretical description of kapellasite, 
which is deep in the Mott phase with well-defined local- 
ized spins, is the Heisenberg Hamiltonian 
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where the exchange integrals J a are defined in Fig. HJa) . 
We use nth order high-temperature (HT) series expan- 
to fit bulk xdc data using J1-J2 or Ji- 
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Jd models, following suggestions that Jd is significant 
in kapellasite [20|. In both cases, Ji comes out ferro- 
magnetic in competition with further-neighbor antiferro- 
magnetic interactions [see red points in Fig.[5][b)], which 
explains the low value of Ocw- Since a nth order HT 
series expansion depends on combinations of up to nth 
points correlation functions that describe the short-range 
arrangements of the spins, it is not surprising that the 
two fits fall in the same phase. The fits are very sta- 
ble, and moving the parameters from the cuboc2 phase 
towards the ferromagnetic one, say, leads rapidly to a 
poor agreement between theory and experiment. The fit 
of the Ji-Jd model (n=ll) is qualitatively better than 
the J1-J2 model and gives a nearest-neighbor interac- 
tion Ji = —15.0(4) K and an "across-hexagon" interac- 
tion Jd = 12.7(3) K. The ferromagnetic nearest-neighbor 
exchange (Ji) highlights the difference between kapel- 
lasite and its polymorph herbertsmithite, and can be 
traced back to the Cu-/i30II~Cu bonding angle being 
~13° smaller in kapellasite [IJ, l21|. The low-T peak in 
the specific heat is also well captured by a HT-series ex- 
pansion using the above Ji-Jd parameters, as shown in 




FIG. 4. (Color online) (a) Wave-vector dependence of the 
inelastic magnetic scattering, S^J^f (Q), integrated over the 
energy range 0.4 < E < 0.8 meV. Symbols are experimental 
data for different temperatures (statistical errors are smaller 
than the symbol size) and the line is the theoretical SBMF 
calculation at T = 5 K. (b) Energy dependence of the imagi- 
nary part of the magnetic dynamic susceptibility, x"{E), in- 
tegrated over the Q range 0.4 < Q < 0.8 A^^ for different 
temperatures. The energy resolution is 0.11 meV. The data 
show deviations from a quasi-elastic Lorentzian (lines) at low 
temperatures. 



Fig.inib). 

Inelastic neutron scattering (INS) measurements were 
performed at 0.1 < T < 300 K on the neutron time-of- 
flight spectrometer INS at the ILL, using incident neu- 
tron energies of 0.82, 3.27, and 20.45 meV. A 93% deuter- 
ated powder sample of mass 3 g was held in an annular 
aluminium can and thermalized by helium exchange gas. 
Standard data reduction 2M gave the neutron scattering 
function S{Q,E), shown in Fig. [TJ The magnetic scat- 
tering, which dominates the nuclear (phonon) scattering 
for the energies and wave vectors studied in this work, is 
quasi-elastic and essentially without dispersion. 

The inelastic magnetic structure factor, 5'™gf (Q), ob- 
tained by integrating S{Q,E) over the energy range 
0.4 < E < 0.8 meV and correcting for the magnetic 
form factor of the Cu^+ ions, is shown in Fig.SJ^a). The 
measured S'-^ff (Q) has initially only a weak tempera- 



incl 

ture dependence, but appears to shift to slightly lower Q 
values for T > 15 K. At even higher temperatures the 
structure gradually disappears, but short-range correla- 
tions persist up to at least 100 K. At low temperatures, 
the broad peak near Q=0.5 A^^ [see Fig.|31Ja)] is centered 
at a much smaller value than in other kagome systems. 
The observed peak position corresponds neither to the T 
point (k = 0) in the Brillouin zone nor to the K point 
k=(l/3,l/3) of the -y/Sx -\/3 structure, but is close to the 
M point, k=(l/2,0). The M point on the kagome lattice 
corresponds to a unique smn arrangement among all pos- 
sible regular spin states [23| of the Hamiltonian ([1]) , called 
cuboc2 (see inset of Fig. [T] and Ref. [2^), and is fully con- 



sistent with the analysis of the HT spin susceptibility and 
the classical phase diagram of the derived Hamiltonian 
[see Fig.[21Jb)]. These two independent measurements, in 
the high- and low-temperature regimes, give strong sup- 
port for that the short-range correlations of kapellasite 
are of non-coplanar cuboc2 type. 

To confirm the relevance of the cuboc2 state for quan- 
tum spins, we have explored the impact of quantum 
fluctuations on the classical system by em ploy ing a 
Schwinger-boson mean-field approach (SBMF) |25[ of the 
Hamiltonian ([T]) , using the values of Ji and J^ extracted 
from xdc and a spin value of S* = 0.5. The powder- 
averaged structure factor S{Q) calculated at T=5 K for 
the same energy range as in the experiment is displayed 
as a line in Fig. IDJa), with only the overall scale fac- 
tor adjusted. The position of the first peak in S{Q) at 
Q = 0.5 A~^ is in good agreement with experiment, as is 
the weak intensity near Q = 1 A^^ . These two features 
are the main signatures of the cuboc2 phase. Classical 
simulations show that replacing Cu by Zn on 27% of the 
kagome sites broaden the structure factor so that the 
peaks at large Qs are partly blurred out while the first 
peak in the powder-averaged S(Q) near Q = 0.5 A~^ ap- 
pears at slightly smaller Qs [26[. Adding quantum fluc- 
tuations to this dilution effect could lead to even better 
agreement with INS data. 

The spin dynamics studied by INS displays a smooth 
continuum of excitations up to at least 2 meV with no 
discernible gap and no sign of damped spin waves, see 
Fig. [T] We extract the energy dependence by integrat- 
ing S{Q, E) over the wave vector range 0.4 < Q < 0.8 
A~^, where the magnetic scattering is strong and the nu- 
clear scattering negligible, and divide by the temperature 
factor 1 — exp(—E/kBT) to obtain the imaginary part 
of the magnetic dynamic susceptibility, x"(£'), shown in 
Fig. Ub). At high temperatures, x"{E) is well described 
by a quasi-elastic Lorentzian, x"{E) = x'ET/[E'^ + F^), 
similar to the classical {S = 5/2) KAFM deuteronium 
jarosite [27|. For T > 5 K, the line width (inverse re- 



laxation rate) increases as T{T) ex T^/^ with tempera- 
ture while the staggered static susceptibility decreases 
as x'{T) ex T~^/^. At lower temperatures, below 5 K, 
the shape of x"{E) deviates from a single Lorentzian, 
which signals the onset of quantum spin-liquid correla- 
tions. However, kapellasite does not show the i?~" be- 
havior of pyrochlore slabs [28| nor the E/T scaling of 
herbertsmithite [J, |29[. This gives further support for 
kapellasite not being close to a quantum critical point. A 
full understanding of the spin liquid dynamics of kapel- 
lasite is not yet at hand and will be discussed further 
below. 

Further insight into the low-T state of kapellasite is 
provided by our /xSR experiments, which were performed 
on a non-oriented powder sample in zero field (ZF) and 
longitudinal applied field (LF) at the ISIS and PSI facili- 
ties down to 20 niK. The evolution of the muon polariza- 
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FIG. 5. (Color online) (a) Time dependent muon polarization 
in zero field (open symbols) and longitudinal field (closed sym- 
bols) demonstrate the absence of spin freezing. Solid lines are 
fits to Eq. ^ while the dashed line represents the expected 
LF polarization at 100 G for frozen spins at 20 mK. (b) Tem- 
perature dependence of the relaxation rate 1/Ti from /^SR 
measurements (circles) and SBMF calculations (red line). 



tion is plotted in Fig. [5ja), both in zero field and under 
iJ = 100 G from 4 K to 20 mK. The relaxation of the 
muon polarization can be fitted in the whole T range to 
the sum of the local responses of two muon sites, known 
to be located nearby OH^ and Cl^ 3, 



Pit) = IPou{t)e 



'AoH* 



(1 - f)Kci{t)e- 



(2) 



where / ^ 0.86 is the T-independent fraction of muons on 
the 0H~ site and Aqh/ci is the dynamical electronic re- 
laxation rate probed at the 0H~/C1^ site. Both Poii{t) 
and the Kubo-Toyabe function Kc\{t) are related to the 
nuclear static fields [3| which are T independent and were 
fixed to their T = 4 K values. 

Any spin freezing can be ruled out based on the follow- 
ing arguments. First, in a fully static scenario, the low-T 
relaxation should be ascribed to a distribution of static 
internal fields with a width /S.H ^ 6 G which should be 
completely decoupled in an applied LF of 100 G. The 
dashed line in Fig. [5ja) shows that this is clearly not the 
case. Second, the ZF polarization at 20 mK relaxes to 
zero, not showing the long-time 1/3-component charac- 
teristic of frozen moments. Finally, the T dependence 
of the ZF polarization is well described by Eq. ^ in 
the whole temperature range from 20 K to 20 mK with 
only one adjustable parameter, XqHi to account for the 
Cu spin dynamics. No electronic relaxation (e"'^'^'* ^ 1) 
was found at the CI site, which is further away from the 
magnetic planes. 

To study the slow spin dynamics, an applied field of 
100 G along the /i+ polarization decouples the static nu- 
clear magnetism so that Eq. ([2]) becomes P{t) = /e^^* + 
(1 — /), and the spin-lattice relaxation rate 1/Ti = Xqh 
is probed directly [Fig. Elb)]. Above 4 K, the 5" = 1/2 



kagomc spins are in a rapidly fluctuating regime. From 
3 K to 0.8 K we observe a slowing down of spin fluctua- 
tions, which occurs in the same T-range as the change in 
dynamics in neutron experiments and coincides with the 
low-T specific heat peak. This indicates an enhancement 
of the short-range correlations at low T, and the corre- 
sponding energy scale is quite low as a result from a com- 
petition between the exchange interactions. Finally, per- 
sistent fluctuations are observed down to 20 mK. This be- 
havior points to a gapless spin-liquid state with an upper 
bound of a hypothetical spin gap of 20 mK, i.e., 10"'^ Ji. 
Within our SBMF approach, the /xSR spin-lattice relax- 
ation rate can be calculated as 1/Ti ex T[x"{ojl^T)/uji\, 
where ui^ is the muon Larmor frequency. The theoretical 
results, also shown in Fig. [5jb), are in excellent qual- 
itative agreement with the experiment, demonstrating 
that the local dynamics of kapellasite is well captured 
by the SBMF calculations at low frequencies. We note 
that the SBMF calculations for both the muon relaxation 
rate and S{Q^ E) give essentially the same predictions for 
both the Ji-Jd and J1--J2 models. This suggests that the 
importance lies not in the exact values of the exchange 
parameters used, but rather in which part of the phase 
diagram kapellasite is located, and this is unambiguously 
determined by the HT expansion. 

At higher frequencies, our SBMF calculations predict 
well-defined spin waves for T < 3 K with energies around 
0.5 meV. We expect that dilution of the magnetic sites re- 
duces the spin-spin correlation length and broadens the 
modes. However, S{Q, E) would still display reminis- 
cences of spin waves, resulting in a non-monotonous en- 
ergy response for certain Q values (such as the magnetic 
zone boundary), which is clearly not the case in kapel- 
lasite, where the energy dependence of the intensity at 
all Qs decrease monotonously. It thus appears clear that 
our a priori choice of bosonic excitations is too restrictive 
to capture both the continuum of spinon excitations and 
the absence of a gap. This indeed calls for further theo- 
retical explorations. A fermionic chiral approach would 
plausibly lead both to gapless and diffuse excitations and 
to a better description of the inelastic neutron spectrum. 
Zero-temperature fermionic approaches of the pure an- 
tiferromagnetic Heisenberg Hamiltonian [9|, |30|, ISlJ have 
been extremely interesting and similar work should be 
extended to the present model with competing interac- 
tions, with the aim to describe both chiral physics and 
thermodynamic properties. Whether the dilution of Cu 
ions is strong enough to kill the predicted low-T chiral 
phase expected with this local spin arrangement J32i | is 
also an interesting question. 
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